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Local  current  density  is  an  important  parameter  in  battery  modeling,  which  affects  the  performance 
of  lithium-ion  batteries.  In  this  study,  we  take  LiFeP04  cathode  material  as  an  example.  A  simplified 
mathematical  model  has  been  developed  to  study  the  internal  mechanism  of  the  electrode.  According 
to  the  results  of  the  model,  the  local  current  density  distribution  has  a  regular  change  at  different  time 
in  the  discharge  process.  The  parameter  “critical  thickness”  as  an  optimized  variable  has  been  presented 
for  battery  design.  By  qualitative  analysis  to  estimate  the  critical  thickness  under  different  condition,  we 
can  optimize  the  design  parameter  of  a  battery  according  to  the  practical  demand. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  modern  times,  lithium-ion  batteries  are  widely  used  for 
portable  electronic  devices,  EVs  and  HEVs.  Among  various  bat¬ 
tery  materials,  LiFeP04  has  been  considered  as  one  of  the  most 
attractive  cathode  material  for  lithium-ion  batteries  due  to  its  low 
cost,  non-toxicity,  great  thermal  stability  and  flat  charge/discharge 
voltage  plateau.  Padhi  et  al.  [1]  first  reported  LiFeP04  as  a  poten¬ 
tial  cathode  material  for  lithium-ion  batteries  in  1997.  Since  then, 
it  has  been  attracting  considerable  attention  as  a  cathode  mate¬ 
rial  for  lithium-ion  batteries  due  to  its  high  theoretical  capacity 
170mAhg_1,  thermal  stability,  and  low  environmental  impact  [2]. 

However,  these  advantages  are  counteracted  by  the  low  elec¬ 
tronic  conductivity  and  ion  conductivity  of  the  material,  which 
is  considered  to  result  in  ohmic  drops  within  the  electrode.  In 
addition,  its  poor  high-rate  capability  also  limits  the  commercial 
utilizations.  Various  ideas  have  been  proposed  to  solve  these  two 
problems,  including  the  reduction  of  the  particle  size  to  reduce 
transport  length  [3]  and  carbon  coating  on  the  particle  surface  to 
improve  the  conductivity  of  the  solid  phase  [4,5]. 

Battery  modeling  is  an  efficient  tool  in  battery  design  and  test, 
which  is  based  on  basic  physical  and  chemical  principles  [6-8,1 1  ]. 
It  is  the  mathematical  description  of  the  internal  structure  and 
the  charging  and  discharging  behavior,  to  find  the  relationship 
between  the  design  variables  and  the  electrochemical  performance. 
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By  performing  a  specific  mathematical  model,  one  can  predict  the 
electrochemical  performance,  capacity  fade  and  thermal  behavior 
of  a  battery.  Usually,  battery  modeling  is  a  complex  system  coupled 
with  several  nonlinear  differential  and  algebra  equations. 

Doyle  and  Newman  firstly  developed  a  general  model  for 
lithium-ion  batteries  [9],  which  is  based  on  porous  electrode  the¬ 
ory:  the  electrode  is  considered  as  a  compound  which  consists  of 
two  phases:  the  porous  matrices  and  the  electrolyte  within  the 
pores.  The  quantities  at  any  point  across  the  electrode  area  are 
assumed  to  be  continuous  functions  of  time  and  space,  and  the 
electrode  is  treated  as  a  superposition  of  two  continua.  The  model 
is  developed  to  simulate  the  charging  and  discharging  curves  of  bat¬ 
teries,  and  predict  the  parameter  distribution  inside  the  battery,  as 
the  potential  distribution  and  ion  concentration  distribution. 

The  mechanism  of  LiFeP04  electrode  is  considered  to  be  unique, 
because  during  the  charge/discharge  process,  lithium  ions  diffuse 
along  a  one-dimensional  pathway  and  phase  transition.  So  the 
mathematical  modeling  of  a  LiFeP04  cathode  requires  describing 
the  behavior  of  the  phase  transition  between  FeP04  and  LiFeP04. 
Srinivasan  and  Newman  [10]  first  develop  a  model  based  on  the 
shrinking-core  model  presented  by  Padhi.  It  is  coupled  with  two 
particle  sizes  and  fitting  based  on  experimental  data. 

In  this  study,  to  make  a  qualitative  analysis,  we  use  a  simpli¬ 
fied  model  based  on  Doyle’s  research  to  estimate  the  performance 
of  LiFeP04  electrode.  We  ignore  the  phase  transition  process  and 
treat  it  as  solid  diffusion.  In  battery  modeling,  local  current  den¬ 
sity  represents  the  reaction  current  density  on  the  surface  of  the 
particles,  which  is  directly  related  to  the  surface  lithium  ion  diffu¬ 
sion  flux.  According  to  the  results  by  solving  the  model  equations, 
we  notice  that  during  the  discharging  process,  the  local  current 
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Fig.  1.  The  schematic  of  the  battery  modeled  in  this  study,  consisting  of  a  Li-foil 
electrode,  a  porous  LiFeP04  cathode  and  a  separator  between  them.  The  battery  is 
filled  with  electrolyte.  Both  electrodes  are  attached  to  a  current  collector. 


density  distribution  across  the  electrode  area  presents  a  peak  value. 
As  discharging  process  goes  on,  the  peak  moves  from  the  current 
collector  to  the  electrode-separator  interface.  This  phenomenon 
will  be  helpful  to  battery  design. 


2.  Model  developments 


The  schematic  of  the  battery  modeled  in  this  study  is  shown  in 
Fig.  1.  The  battery  consists  of  a  Li  foil  anode  and  a  LiFeP04  cathode 
with  a  porous  separator  between  the  two  electrodes.  The  cathode 
consists  of  a  large  amount  of  solid  spherical  particles.  So  the  porous- 
electrode  theory  can  be  introduced  in  this  model.  We  developed  a 
one-dimensional  model  in  this  study.  The  mathematical  equations 
that  describe  the  processes  are  listed  below. 

In  the  electrode  scale,  according  to  porous  electrode  theory,  the 
current  density  in  the  electrolyte  phase  is  determined  by  a  modified 
Ohm’s  Law: 


h  =  V  • 


— /c2  V  &2  + 


2RTk2 

F 


1  + 


9  In/ 
9  In  c2 


[1  -r+]V(ln  c2) 


(1) 


Flere  the  subscript  *2'  refers  to  the  electrolyte  phase.  @2  is  the 
potential,  c2  is  the  concentration  of  lithium  ion  and  k2  is  the  ion 
conductivity  in  the  electrolyte. 

In  the  solid  phase  (denoted  by  subscript  ‘1'),  the  current  density 
is  determined  by  Ohm’s  law: 


h  =  I-i  2  = 


(2) 


where  I  is  the  applied  current  density.  As  a  consequence  of  the 
assumption  of  electroneutrality  that  the  divergence  of  the  total  cur¬ 
rent  density  is  zero.  For  a  one-dimensional  model,  it  is  expressed 
as 


ii +i2=/  (3) 

According  to  porous-electrode  theory,  the  mass  balance  in  the 
electrode  is 


A  rate  equation  is  used  to  determine  the  dependence  of  the  local 
electrochemical  reaction  rate  on  concentration  and  potential.  The 
Butler-Volmer  equation  is  commonly  used: 

i^,o{exp(S)“exp(i r)}  (5) 

Flere  in  refers  to  the  local  current  density  and  rj  =  -  @2  -  Eocv 

refers  to  the  surface  overpotential.  E0Cv  is  the  open-circuit  poten¬ 
tial  of  LiFePC>4  with  respect  to  a  lithium  reference  electrode,  which 
is  a  function  of  the  state  of  discharge  (SOD),  and  can  be  evaluated 
by  discharging  at  very  low  rate.  Generally  the  state  of  discharge 
can  be  represented  as  an  expression  of  solid  phase  concentration 
Cl :  Ci/c  max  [10].  Cmax  is  the  maximum  concentration  of  lithium 
ions  that  can  be  incorporated  into  the  lattice,  which  is  take  to  be 
22785  mol  m-3  for  LiFeP04.  The  exchange  current  density  i0  is  given 
by: 

i0  =  Fkc°-5c°-5(cmax  -  ci  )0'5  (6) 

where  (cmax  -  C\ )  denotes  the  concentration  of  unoccupied  sites  in 
LixFeP04  and  k  is  the  reaction  rate  constant. 

A  charge  balance  relates  the  divergence  of  the  current  to  the  sur¬ 
face  diffusion  flux  jn,  or  the  local  current  density  in  due  to  reaction: 

ain  =  aFjn  =  V  i2  (7) 

The  boundary  condition  is:  at  the  current  collector  interface 
(x  =  0),  i2  =  0;  at  the  electrode-separator  interface,  i2=l. 

In  the  solid  phase,  the  phase  transition  from  FeP04  (charged)  to 
LiFeP04  (discharged)  during  charge/discharge  is  ignored.  The  pro¬ 
cess  of  lithium  ions  insertion/extraction  is  treated  as  solid  diffusion, 
and  follows  the  second  Fields  Law  in  spherical  coordinates: 

dCs  d2Cs  2D,  dCs  m 

9F=Dl9^  +  ~9F  (8) 

The  boundary  condition  is: 


r  =  0 


i  -  l—  -  -D  — 
Jn  F  Dl  dr  ’ 


r  =  Rs 


The  solid  diffusion  is  based  on  particle  scale,  where  Rs  is  the  par¬ 
ticle  radius.  cs  denotes  to  the  solid  phase  concentration  of  lithium 
ions  in  the  particles.  At  the  surface  of  the  particles,  the  solid  con¬ 
centration  is  corresponding  to  the  solid  phase  concentration  C\  in 
the  electrode  scale.  The  surface  diffusion  flux  of  lithium  ions,  jn,  is 
embedded  into  Eqs.  (4)  and  (5)  in  the  electrode  scale. 

All  equations  ((1),  (2),  (4),  (5),  (7)  and  (8))  are  simultaneously 
solved  using  Newton’s  iterative  method  [12].  The  electrode  domain 
is  discretized  using  finite  difference  method.  Six  unknowns  are 
involved  in  this  model,  including  the  concentration  in  the  solid 
and  electrolyte  phase  C\  and  c2,  the  potential  in  the  two  phases 
and  02,  the  current  density  in  the  electrolyte  phase  z2,  and  the 
local  current  density  in,  respectively.  At  each  time  step,  the  iter¬ 
ation  commences  by  providing  initial  guesses  for  all  unknowns, 
until  the  convergence  is  reached.  Time  stepping  moves  on  until  the 
cut-off  voltage  is  reached.  From  the  results  we  can  obtain  the  dis¬ 
tribution  of  all  the  unknowns  across  the  electrode  area  at  different 
time.  The  parameters  used  in  the  model  are  listed  in  Table  1  and  the 
open-circuit  potential  curve  is  shown  in  Fig.  2,  which  is  obtained 
by  discharging  the  cell  at  0.02  C. 


3.  Results  and  discussions 


£ 


dc2 

dt 


=  V.(DVc2)+^(  1-4) 


(4) 


During  the  discharge  process,  lithium  ions  migrate  and  diffuse 
from  anode  to  cathode  across  the  separator  and  insert  into  the 
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Table  1 

Parameters  used  in  the  model. 


Parameter 

Value 

Diffusion  coefficient  in  LiFeP04,  D 1 

9  x  10-18  m2  s_1 

Solid  phase  conductivity,  K\ 

5  x  10~3  Snrr1 

Electrolyte  conductivity,  K2 

[10] 

Particle  radius 

4  pan 

Initial  salt  concentration 

1M 

Reaction  rate  constant,  k 

2  x  10-11 

Separator  thickness 

25  pan 

Transfer  coefficient,  a 

0.5 

Ion  transfer  number  for  lithium  ions,  t+ 

0.39 

active  material  particles.  The  electrode  reaction  occurs  on  the  sur¬ 
face  of  the  particles,  resulting  in  a  surface  diffusion  flux  of  lithium 
ions  jn  into  the  particles.  Based  on  the  Butler-Volmer  equation 
(Eq.  (5)),  the  local  current  density  in  is  produced  by  the  electrode 
reaction,  and  in  the  meantime  it  is  described  as  a  flux  boundary 
condition  in  the  solid  diffusion  process  (Eq.  (8)). 

In  the  solid  diffusion,  what  we  care  about  is  the  relation  between 
the  surface  diffusion  flux  jn /local  current  density  in  and  the  sur¬ 
face  concentration  of  lithium  ions  c\.  By  discretizing  Eq.  (8)  using 
backward  differences  in  the  radial  coordinates,  a  following  form  of 
matrix  equation  can  be  obtained: 

Ac*  =  x  C5-1  +  bk  (9) 

where  A  is  a  tridiagonal  coefficient  matrix,  cs  is  the  concentration 
vector  inside  the  particles,  b  is  a  vector  relative  to  the  local  current 
density  in  and  k  denotes  the  present  time  step  while  k  -  1  denotes 
the  previous.  By  solving  this  equation  using  Thomas  algorithm,  an 
expression  related  local  current  density  to  surface  concentration 
on  the  particles  can  be  obtained: 

4  =  pcj-1  +  qi'n  (10) 

where  p  and  q  are  coefficients,  which  vary  at  different  time  step 
and  different  position  in  the  electrode  domain.  In  the  discharge 
process,  as  the  electrode  reaction  produces  surface  flux,  lithium 
ions  insert  into  the  particles  continuously.  As  a  result,  the  surface 
concentration  of  lithium  ions  increases  and  the  state  of  discharge 
also  rises.  While  the  surface  concentration  rises  up  to  a  value  closed 
to  cmax,  namely  the  local  SOD  is  close  to  1,  the  lithium  ions  will  be 
difficult  to  be  inserted. 


V 


We  simulate  the  local  current  density,  surface  concentration  and 
local  state  of  discharge  distribution  across  the  cathode  area  at  1  C 
discharge.  To  make  a  general  discussion,  the  cathode  thickness  is 
set  to  be  100  fxm.  The  curves  at  different  time  in  the  discharge  are 
shown  in  Figs.  3-5.  In  Fig.  3,  a  peak  exists  in  the  local  current  density 
curve  which  moves  from  current  collector  to  the  cathode-separator 
interface  as  discharge  proceeds.  At  the  beginning  of  discharge,  the 
peak  presents  at  the  current  collector  interface,  namely  the  local 
current  density  at  the  boundary  is  much  higher  than  that  in  the  inte¬ 
rior,  indicating  a  much  faster  reaction  rate  near  the  current  collector 
region.  Therefore,  the  surface  concentration  near  the  current  collec¬ 
tor  increases  very  fast  as  shown  in  Fig.  4.  As  the  reaction  proceeds, 
the  surface  concentration  in  this  region  will  rise  up  to  be  closed 
to  maximum  concentration  and  the  particles  will  approximately 
reach  the  fully  discharged  state  in  a  short  time  (Fig.  5). 

The  battery  is  discharged  under  galvanostatic  conditions.  By 
charge  balance,  the  amount  of  lithium  ions  inserted  into  particles 
will  be  constant  per  unit  time,  so  the  integral  of  local  current  den¬ 
sity  across  the  cathode  area  is  also  a  constant  relative  to  the  applied 
current  density  (Eq.  (7)).  When  the  particles  close  to  the  current 
collector  approximately  reach  their  fully  discharged  state,  lithium 
ion  can  hardly  continue  to  be  inserted,  resulting  in  a  very  low  local 
current  density  in  this  region  during  the  later  discharge  process.  To 
satisfy  the  charge  balance  condition,  the  peak  moves  to  the  leftward 
region  until  the  local  state  of  discharge  approximately  equals  to  1 
(as  shown  in  Fig.  3).  This  process  repeats  until  the  end  of  discharge 
is  reached. 

In  Figs.  4  and  5,  the  surface  concentration  and  state  of  discharge 
distributions  at  different  time  represent  the  same  trend.  The  parti¬ 
cles  close  to  the  current  collector  boundary  are  firstly  approaching 
the  fully  discharge  state.  As  a  consequence,  we  can  conclude  that 
high  local  current  density  means  high  diffusion  flux  and  reaction 
rate  which  results  in  fast  growth  rate  of  surface  concentration  and 
state  of  discharge. 

As  the  curve  at  t  =  3575  s  shown  in  Fig.  3,  we  can  find  that  at 
the  end  of  discharge,  the  peak  of  the  local  current  density  has 
not  reached  the  separator-cathode  boundary.  Correspondingly,  the 
surface  concentration  in  the  region  close  to  separator  also  has  not 
reached  the  maximum  concentration  and  the  local  state  of  dis¬ 
charge  is  far  less  than  1,  which  indicates  that  the  active  material 
has  not  been  adequately  utilized.  That  is  because  the  stop  condition 
is  determined  by  cut-off  voltage  in  both  simulation  and  practice. 
When  cut-off  voltage  is  reached,  the  utilization  of  active  material 
may  not  be  1 00%.  As  a  result,  to  optimize  battery  design  and  achieve 
better  performance,  it  is  essential  to  choose  appropriate  value  for 
some  design  parameters. 

In  electrode  design,  porosity  e,  electrode  thickness  h,  compacted 
density  pc  and  surface  density  ps  (the  mass  of  active  material  per 
unit  electrode  area)  are  four  important  design  parameters.  For  a 
specific  electrode  material  (the  intrinsic  density  has  been  fixed),  if 
the  component  of  the  slurry  is  given  (including  the  active  mate¬ 
rial,  conductor  and  polymer  binder),  the  four  parameters  can  be 
approximately  related  by  the  following  equations: 


3.2 


£  —  \  £act  Scond  & bind 


Shind  —  1  Pc  f 


■  + 


xcond 


\  Pact  xactPcond 
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xbind 


ind  \ 
°bind  J 


X-actPbind  > 

(11) 
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Fig.  2.  The  equilibrium  potential  curve  for  LiFeP04.  It  is  obtained  by  discharging  at 
0.02  C. 


where  pact,  pconc/,  Pbind  are  the  intrinsic  density  of  active  material, 
conductor  and  binder,  xacu  xcond .  xbind  denote  the  mass  fraction  of 
the  three,  and  sact,  £cond,  £bind  denote  the  volume  fraction  of  them  in 
the  electrode,  respectively.  If  the  porosity  and  the  electrode  thick¬ 
ness  are  fixed,  the  other  two  parameters  can  also  be  determined 
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Distance  across  cathode  x  (Jim) 


Fig.  3.  The  local  current  density  distributions  at  different  time  during  the  discharge  at  1  C  rate.  The  cathode  thickness  is  set  to  be  100  |jim.  The  abscissa  denotes  the  cathode 
thickness:  the  left  side  represents  the  separator-electrode  interface  and  the  right  side  the  current  collector  boundary.  The  numbers  in  the  legend  denote  the  time  of  discharge 
in  seconds.  The  value  of  local  current  density  is  negative  which  means  the  charges  flow  into  the  particles. 


from  the  equations  above.  As  a  consequence,  we  need  to  analyze  the 
effect  of  two  parameters:  porosity  and  electrode  thickness,  which 
are  also  introduced  as  parameters  in  our  model,  on  the  performance 
of  a  battery. 

To  achieve  the  optimum  utilization  of  the  active  material, 
we  should  select  an  appropriate  electrode  thickness.  From  the 


discussion  above,  we  may  derive  a  ‘critical  thickness’  as  the  opti¬ 
mum  thickness.  It  is  determined  by  the  location  of  the  peak  in 
the  local  current  density  distribution.  When  the  cut-off  voltage  is 
reached,  if  the  peak  value  of  local  current  density  just  arrives  at 
the  electrode-separator  interface,  namely,  the  active  material  in 
the  whole  electrode  has  been  adequately  utilized,  this  thickness 


Fig.  4.  The  surface  concentration  distributions  at  different  time  during  the  discharge  at  1  C  rate. 
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Fig.  5.  The  state  of  discharge  distributions  at  different  time  during  the  discharge  at  1  C  rate.  The  fully  discharged  state  moves  from  the  current  collector  to  the  interior  as 
discharge  proceeds. 


is  considered  to  be  the  “critical  thickness”  under  the  condition.  If 
this  electrode  thickness  is  greater  than  critical  thickness,  when  the 
cut-off  voltage  has  been  reached,  the  peak  has  not  arrive  at  the 
electrode-separator  interface  yet,  indicating  there  are  some  redun¬ 
dant  active  material  which  has  not  been  utilized.  Otherwise,  if  the 
thickness  is  less  than  critical  thickness,  before  the  cut-off  voltage 
is  reached,  the  peak  has  already  reached  the  boundary,  indicating 
that  the  active  material  may  be  inadequate. 


In  other  words,  to  calculate  the  critical  thickness,  two  stop  con¬ 
ditions  should  be  satisfied  at  the  same  time:  the  cut-off  voltage 
is  reached  and  the  peak  of  local  current  density  distribution  just 
arrives  at  the  separator-cathode  boundary.  For  a  specific  electrode 
component,  the  critical  thickness  is  related  to  the  porosity  of  the 
electrode  and  the  discharge  current  density/discharge  rate. 

In  the  following  discussion,  we  will  make  a  qualitative  analy¬ 
sis  about  the  effect  of  porosity  and  discharge  rate  on  the  critical 
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Fig.  6.  The  local  current  density  distributions  at  different  time  during  the  discharge  at  0.5  C  rate.  The  cathode  thickness  and  the  porosity  are  set  to  be  0.45  and  80  |jim, 
respectively.  The  peak  approximately  moves  to  the  separator  when  cut-off  voltage  is  reached,  indicating  the  critical  thickness  is  80  |jim  under  this  condition. 
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Fig.  7.  The  surface  concentration  distributions  at  different  time  during  the  discharge  at  0.5  C  rate.  The  porosity  and  thickness  are  0.45  and  80  |jim.  When  the  thickness  is 
designed  to  be  the  critical  thickness,  the  active  materials  are  adequately  utilized  at  the  end  of  discharge. 


thickness.  We  take  three  porosity  values  (0.45,  0.4,  0.35)  and  four 
discharge  rates  (0.5  C,  1  C,  2  C,  5  C),  and  calculate  the  local  current 
density  distribution  at  different  time  of  various  electrode  thickness, 
in  order  to  find  the  critical  thickness. 

When  the  porosity  equals  to  0.45,  the  local  current  density 
and  surface  concentration  distributions  at  0.5  C  are  shown  below 
(Figs.  6  and  7).  From  the  figures  we  can  find  that  the  critical 
thickness  is  80  p,m  at  this  discharge  rate.  When  cut-off  voltage  is 
reached,  the  peak  of  the  local  current  density  curve  approximately 
reaches  the  boundary  at  the  same  time.  The  calculation  process 
is  repeated  then  the  critical  thickness  under  various  conditions  is 
obtained. 

All  the  results  are  summarized  in  Table  2.  We  can  infer  if  the 
porosity  has  been  fixed,  as  the  discharge  rate  grows,  the  critical 
thickness  becomes  less;  if  the  discharge  rate  is  fixed,  as  the  porosity 
reduces,  the  critical  thickness  also  decrease  (Table  2). 

It  is  only  a  qualitative,  not  a  quantitative  analysis,  but  still  has 
great  significance  on  battery  design.  There  is  an  experienced  rule 
in  the  battery  industry:  to  make  better  battery  for  higher  discharge 
rates,  the  electrode  thickness  should  be  smaller.  Our  results  provide 
some  theoretical  support  on  this  rule.  Once  demand  of  design  has 
been  confirmed,  we  can  manufacture  the  corresponding  batteries 
with  the  optimized  thickness  which  matches  compacted  density 
(porosity)  and  discharge  rates. 


Table  2 

Critical  thickness  (|jim)  for  LiFeP04  of  different  porosity  and  discharge  rates. 


Porosity 

Discharge  rates  (C) 

0.5 

1 

2 

5 

0.45 

80 

70 

55 

30 

0.4 

75 

60 

50 

20 

0.35 

70 

55 

40 

15 

4.  Conclusions 

Local  current  density  is  an  important  parameter  in  battery 
modeling,  which  represents  the  surface  ion  flux  and  relates  to 
the  electrochemical  reaction  rates.  Based  on  the  solution  of  the 
mathematical  model  we  can  estimate  the  local  current  density  dis¬ 
tribution  across  the  electrode  area.  At  different  time  during  the 
discharge  process,  the  curve  appears  the  different  trend.  A  peak  is 
discovered  on  the  curves  and  it  moves  on.  By  qualitative  analysis  we 
get  a  parameter  called  “critical  thickness”  which  may  be  important 
for  battery  design.  In  further  research  we  will  make  quantitative 
analysis  to  get  accurate  results. 
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Appendix. 


List  of  symbols 

a  specific  interfacial  area  per  unit  volume  of  electrode 

(m-1) 

c2  salt  concentration  in  the  electrolyte  (mol  m-3 ) 

Ci  concentration  of  lithium  in  solid  phase  (mol  m-3 ) 

Cmax  maximum  concentration  of  lithium  in  solid  phase 

(molnrr3) 

Di  diffusion  coefficient  of  lithium  ion  in  solid  phase  (m2  s-1 ) 
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D?  diffusion  coefficient  of  lithium  ion  in  the  electrolyte 

(m2  s-1) 

F  Faraday’s  constant  (96,487  C  mol-1 ) 

h  electrode  thickness  (m) 

I  applied  current  density  (A) 

i0  exchange-current  density,  referred  to  the  active  surface 

area  (Am-2) 

in  local  current  density  at  the  particle  surface,  referred  to 

the  active  surface  area  (Am-2) 

Rs  particle  radius  (m) 

t+  ion  transport  number  of  the  electrolyte 

R  universal  gas  constant  (8.314 J  mol-1  K-1 ) 

£  porosity,  the  volume  fraction  of  electrolyte  in  the  elec¬ 

trode 

0i  potential  of  the  solid  phase  (V) 

02  potential  of  the  liquid  phase  (V) 

rj  surface  overpotential  (V) 

a  transfer  coefficient 

pc  compacted  density,  mass  of  active  materials  per  unit  vol¬ 
ume  of  the  electrode  (kgm-3) 


ps  surface  density,  mass  of  active  materials  per  unit  area  of 
the  electrode  (kgm-3) 
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